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O . ABSTRACT 

o 



We report on an XMM-Newton Reflection Grating Spectrometer observation of the 
^S| ■ black hole candidate and Galactic microquasar GRS 1758—258. The source entered a 

peculiar "off /soft" state in late February, 2001, in which the spectrum softened while 
the X-ray flux - and the inferred mass accretion rate - steadily decreased. We find 
j;^ I no clear evidence for emission or absorption lines in the dispersed spectra, indicating 

C^ ' that most of the observed soft flux is likely from an accretion disk and not from a cool 

plasma. The accretion disk strongly dominates the spectrum in this lower-luminosity 
state, and is only mildly recessed from the marginally stable orbit. These findings may 



o 

:^ 

o . be difficult to explain in terms of advection-dominated accretion flow, or "ADAF" 

models. We discuss these results within the context of ADAF models, simultaneous 
two-flow models, and observed correlations between hard X-ray flux and jet production. 



k> 1. Introduction 

. P_ , The source GRS 1758-258 was discovered with the GRANAT /SIGMA hard-X-ray/soft-7-ray 

telescope (Sunyaev et al. 1991). It is often referred to as a "twin" source with IE 1740.7—2942. 
Both are in the vicinity of the Galactic Center region, are strong emitters at hard-X-ray/soft-7-ray 
energies, and display relativistic jets in the radio band (Mirabel et al. 1992; Rodriguez, Mirabel, 
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& Marti 1992; Mirabel 1994). The hard spectrum of GRS 1758-258 extends to 300 keV and is 
similar to that of Cygnus X-1, and on that basis it is considered a black hole candidate (for a 
review of BHCs, see Tanaka & Lewin 1995). 

The soft X-ray component in this source is usually very weak when the hard X-ray flux is 
at a typical strength. The soft component was too weak to be positively detected with ASCA 
(Mereghetti et al. 1997). However, a soft component below 2 keV, modeled with a power-law 
(FpL ~ 3) was observed with ROSAT in 1993 when the hard power-law flux was less intense 
(Mereghetti et al. 1994). Similarly, the soft X-ray component was detected with XMM-Newton 
during a short, softer episode in September, 2000 (Goldwurm et al. 2001). In this observation, the 
soft component could be fit with a simple blackbody spectrum (kT = 0.32 keV). 

In late February of 2001, observations with RXTE revealed a sharp drop in the hard 
(>10 keV) flux from GRS 1758—258, with no corresponding drop in the soft flux (Smith et 
al. 2001a, b,c). Smith et al. (2001c; hereafter, SHMS) report a steady decline in the 3-25 keV 
flux for 50 days following 27 February 2001. The soft component is clearly detected with the 
RXTE /PCA, and can be fit by both simple blackbody and multicolor disk (MCD; Mitsuda et 
al. 1984) blackbody models. This state is qualitatively similar to the "off" state observed by 
Gi^^Af^T/SIGMA between fall 1991 and spring 1992 (Gilfanov et al. 1993), during which the 
hard X-ray flux of GRS 1758—258 decayed beneath detection limits. We therefore refer to the 
present X-ray spectral state as the "off/soft" state. 

We requested a XMM-Newton target of opportunity (TOO) observation of GRS 1758—258 for 
three principal reasons: to better understand the nature of the soft component, to understand the 
implications of the off /soft state for current accretion flow models, and to place the off /soft state 
within the context of emerging connections between spectral flux components and jets in BHCs 
(Fender 2001). Herein we report the results obtained with the XMM-Newton Reflection Grating 
Spectrometer (RGS). 



2. Observations and Data Analysis 

Following our TOO request, GRS 1758-258 was observed for ~22 ksec starting at UT 09:48:12 
on 22 March, 2001. RGSl and RGS2 were operated in "spectral-^Q" mode. The EPIC-MOS 
cameras were operated in "partial window, 100x100" mode, and the EPIC-PN camera in "large 
partial window, 200x384" mode. For the OM, "grism 1," optimized for the ultraviolet band, was 
used. 

Accurate analysis of bright sources like GRS 1758—258 requires a careful treatment of pile-up 
when dealing with the EPIC data. This in turn requires a good knowledge of the mirror response; 
the EPIC calibration team is developing a description of the in-flight calibration. The results 
presented in this work are therefore restricted to analysis of the RGS data, and results from the 
EPIC data analysis will be presented in future work. At the time of writing, the only previous 



X-ray binary observations reported on with XMM-Newton are EXO 0748—676 (Bonnet-Bidaud et 
al. 2001; Cottam et al. 2001) and GRS 1758-258 (observation in September 2000, Goldwurm 
et al. 2001). Instrumental performance and systematic effects will become more clear as more 
sources with a strong soft continuum flux are observed. 

We have used the spectral data and background files produced by the standard pipeline 
processing system at the XMM-Newton Survey Science Center (SSC) using Science Analysis 
Software (XMM-SAS) version 5.0. We estimate that fewer than 1% of the counts in the RGS 
spectra are from the nearby source GX 5—1. 

To obtain improved statistical constraints on fitted models and to ensure the accuracy of 
Poisson statistics, we rebinned the background-subtracted first-order spectra from RGSl and 
RGS2 by requiring a minimum of 20 counts per bin. After rebinning, between 0.35-0.60 keV the 
flux bin errors include zero. The flux bins steadily increase and the errors are inconsistent with 
zero above 0.60 keV. We therefore adopt 0.60 keV as the lower limit of our fitting range. Similarly, 
above 2.3 keV the flux bins have large errors and are consistent with zero, so we adopt 2.3 keV as 
the upper limit to our fitting range (the full energy range of the RGS is 0.35-2.50 keV). 

Using XSPEC version 11.0 (Arnaud & Dorman, 2000), the RGSl and RGS2 spectra are 
fitted simultaneously to account for any cross-calibration uncertainties. An overall normalization 
constant is allowed to float between the data obtained from each. Values obtained for this constant 
are approximately 0.95 for all fits. 



3. Results 

The results of our joint fits to the RGS spectra are detailed fully in Table 1. In each of the 
six models, the continuum components are multiplied by a model for photoelectric absorption 
("wabs," in XSPEC). Model 6 is a single-component model consisting of a simple power-law. 
Although statistically only marginally worse than models dominated by a thermal component 
k)d 1^ = 1.778,1^ = 1251), the power-law index is abnormally high (TpL = 5.2 it 0.2), and 
inconsistent with values obtained for other BHCs regardless of spectral state (for a review, see 
Tanaka & Lewin 1995). We therefore conclude that models dominated by a thermal component 
are required to fit the data from GRS 1758—258 in the off/soft state, and concentrate on these. 

In models which include thermal and power-law components, the power-law index could not 
be constrained. In order to be consistent with models widely fit to BHCs, and especially those 
recently fit to GRS 1758—258, we fix the power-law indices in models 1-5 to those measured by 
SHMS via fits to RXTE /VCA data obtained on 12-13 March, 2001 in the 3-25 keV band using 
the same models. The fit results for models 1-5 indicate that the power-law component is not 
required in this energy range as a zero normalization is preferred. We quote 90% upper-limits 
for the 0.6-2.3 keV contribution from a power-law component in Table 1. In all cases, the 
power-law contributes less than 2% of the total flux in this band. This is consistent with Chandra 



measurements in the 1-10 keV band during the off /soft state (Heindl and Smith 2001). 

Models 1 and 3 consist of MCD and simple blackbody components, respectively, in 
combination with a power-law component {x^ /v = 1.621, x^ /u = 1.587, respectively, v = 1252). 
The data/model ratios obtained for models 1 and 3 are nearly identical (see Figure 1), and neither 
indicates the presence of any significant line emission (indeed, line emission is not indicated 
regardless of continuum model). In part, the poor x^ statistics obtained for these models (indeed, 
for all models) may be due to the approximately 5% normalization discrepancy between RGSl 
and RGS2. A second contributing factor may be the relatively smaller RGS effective area 
above ~2 keV; we included this region after rebinning to better-constrain any power-law flux. 
When the characteristics of the RGS are more clearly known for strong continuum sources like 
GRS 1758—258, it may be possible to improve on the fits. 

Models 2 and 4 are duplicates of models 1 and 3 respectively, except in models 2 and 4 we 
have fixed the neutral hydrogen column density (Njy) to 1.5x10^^ atoms/cm^ as per fits to ASCA 
data of GRS 1758-258 reported by Mereghetti et al. (1997). Model 5 is like 3 and 4 in that 
simple blackbody and power-law components are used to model the continuum, but in this model 
Nh is fixed to 1.74x10^^ atoms/cm^, as measured via fits to XMM-iVeu;ton/EPIC-MOS data on 
GRS 1758—258 obtained in September, 2000 (Goldwurm et al. 2001). There are no reasons to 
doubt the value of N// we measure due to instrumental effects (e. g., GCD pile- up, as the RGS is 
a dispersive spectrometer). Rather, these models are fit to allow for comparisons to the results of 
SHMS and Goldwurm et al. (2001) as directly as possible. 

Via the MCD model, SHMS obtain a color temperature of kT = 0.464 ± 0.007 keV. Via 
model 1 (Hh free), we find a color temperature of kT = 0.34 it 0.01 keV. In both measurements, 
the errors are 90% confidence errors; these color temperatures are significantly different. Via 
model 2 {Hh fixed at SHMS value), we find a color temperature of kT = 0.60 it 0.01 keV. 
Again, this is significantly different from the value measured by SHMS. However, the fit 
with model 2 {x^ /f = 2.227, u = 1251) is statistically worse than that obtained via model 1 
{x^/u = 1.621,1/ = 1252). 

SHMS report a simple blackbody temperature of kT = 0.395 ± 0.006 keV. Fits with model 
3 (N// free) indicate kT = 0.286 ± 0.007 keV, and via model 4 {^h fixed at SHMS value) 
kT = 0.378l']Q'gQ4 keV. The temperature measured via model 4 is only slightly lower than that 
measured by SHMS, and only marginally inconsistent at 90% confidence. However, the fit obtained 
with model 3 is statistically preferred and the temperature is significantly lower. 

The blackbody temperature obtained via fitting with model 5, wherein N// is fixed to the value 
measured by Goldwurm et al. (2001), is kT = 0.332l[]:[![j^ keV. They measure kT = 0.32 ±0.02 keV 
- these temperatures are consistent at 90% confidence. Statistically, model 3 is a marginally better 
fit, and for this model the measured temperature is significantly lower. 

We attempted to describe the spectrum in terms of a diffuse plasma by fitting with the 
"mekal" and "raymond" models within XSPEC. Fits which are statistically only marginally worse 



than those with models 1-6 are obtained, but only if the elemental abundances are allowed to 
assume values less than 0.2% of solar values. As this physical scenario is extremely unlikely, we do 
not further consider these models. 

To test for the presence of narrow emission or absorption lines, we apply the continuum 
we measured via model 1 in fits to the rebinned data to the RGSl and RGS2 spectra at full 
instrumental resolution. We selected model 1 instead of model 3 as we regard the MCD model to 
be more physical for BHCs than a simple blackbody model (the difference in the continuum shape 
from model three should be very small). In Figure 2, we plot the 90% confidence upper-limits 
on the strengths of emission or absorption features with widths less than or equivalent to the 
first-order RGS resolution. For visual clarity, we have rebinned the data in this plot by a factor of 
10. It is clear from this plot that any narrow line features in the 0.6-2.3 keV range are very weak 
(equivalent width ~l-5 eV for most of the bandpass) when GRS 1758—258 is in the off/soft state. 

Finally, we examine the lightcurve obtained during this observation for the full RGS energy 
range (see Figure 3). Standard Shakura-Sunyaev (1973) accretion disks should vary on the 
viscous timescale through the disk (weeks). In sources wherein a power-law is observed (and a 
coronal volume therefore implied), significant short-timescale flux variations (seconds-hours) are 
observed (see Tanaka & Lewin 1995). Such rapid variations - both aperiodic and quasi-periodic - 
can be described in terms of magnetic flaring between the disk and corona or within the corona 
itself (see, e.g., Di Matteo, Celotti, k Fabian 1999; Merloni, Di Matteo, & Fabian 2000). The 
lightcurve shown in Figure 3 is extremely steady, supporting an accretion disk interpretation for 
the spectrum we have observed. 



4. Discussion 

Although we have observed the soft X-ray component in the off /soft state of GRS 1758—258, 
no clear emission or absorption features have been observed in the joint spectra of RGSl and 
RGS2 (see Figures 1 and 2). Observations of BHCs IE 1740.7-2942 (Cui et al. 2001) and 
SS 433 (Marshall et al. 2001) with Chandra have found evidence for extended regions around 
the X-ray sources which are too large to be an accretion disk but too small to be a supernova 
remnant. A central goal of our observation was to examine the nature of such an extended region 
in GRS 1758—258, if one exists. Since optically-thin, cooling plasmas produce strong line features, 
it is unlikely that the soft component we have measured is due to such a region. The improbable 
abundances required by flts with models for diffuse plasmas support this interpretation. 

The jets from GRS 1758—258 observed in the radio band might be another potential source for 
X-ray line emission. Tentative evidence for lines from a jet are found in the Chandra observation 
of IE 1740.7—2942 noted above; lines are clearly detected in the spectra of SS 433. The absence 
of lines that might be attributed to a jet in GRS 1758—258 is consistent with a picture in which 
jets are extinguished in the soft state in BHCs (Fender 2001). Recently, it has also been proposed 



that jets might produce the power-law component observed in BHCs (MarkofF, Falcke, & Fender 
2001). The observed absence of hnes consistent with a jet, coupled with the diminished strength 
of the power-law component in the off/soft state, may support this model. Stronger conclusions 
must await the publication of radio data gathered during the off/soft state of GRS 1758—258. 

The soft X-ray component we have measured is likely due primarily (but not necessarily 
entirely) to an accretion disk. This interpretation is bolstered by the lack of significant variability 
in the lightcurve obtained for this observation (see Figure 3). We measure a disk temperature 
which is broadly consistent with results obtained with the RXTE /FCA by SHMS on 12-13 March 
2001, and with a short soft state only five months previous to our observation (Goldwurm et al. 
2001). From Table 1, it is clear that the temperature of the soft component and the 0.6-2.3 keV 
flux are dependent upon the value of Nj:/ which is used. Based on the results of SHMS, and our 
finding that a power-law component contributes less than 2% of the flux in the 0.6-2.3 keV range, 
it is possible that in the extended off/soft state the hard power-law flux merely turns-off and 
reveals the accretion disk. 

The off /soft state is indicated by a steady decrease in the 3-25 keV flux from GRS 1758—258, 
and by a far softer spectrum than is observed when the source is in the more typical low/hard 
state (SHMS 2001). Such a decrease in flux is usually interpreted in terms of a decreasing rh. 
Advection-dominated accretion flow (ADAF) models adapted to BHCs (Esin, McClintock, & 
Narayan 1997; Esin et al. 1998) predict that a decrease in the mass accretion rate (rh) should be 
accompanied by spectral hardening. The spectral hardening may occur as the inner accretion disk 
is replaced by a hot, optically-thin region which may Comptonize seed photons and produce a 
power-law spectral component. In sharp contrast, the spectrum of GRS 1758—258 softens as the 
inferred rh decreases. 

If the distance to a source and its inclination are known, the MCD model provides a measure 
of the temperature and inner radial extent of the accretion disk. In fact, the derived temperature 
and radius may be distorted by the Comptonization of disk photons by a coronal volume or a disk 
atmosphere. Shimura and Takahara (1995) have reported that a simple correction factor f^oi can 
account for this effect; the effective disk temperature is given by Tg// = Tcoi/fcoi, ^^^ the effective 
radius is given by v^ff = fcoi'^T^coi- For BHCs, a typical value is fcoi = 1-7 (see, e. g., Sobczak et al. 
1999). 

In Figure 4, we present constraints on the inner radial extent of the accretion disk. The 
best-fit color temperature and radius measured with MCD models in Table 1 are used. We assume 
a distance of 8.5 kpc to GRS 1758—258 based on its central Galactic position and column density. 
The inclination of GRS 1758—258 has not yet been measured, however the detection of jets in the 
radio band makes it unlikely that the system is seen face-on. We therefore assume an intermediate 
value for the inclination angle {9ind = 45). We plot constraints on the inner radii calculated using 
no color correction {fcoi = 1-0) and the typical value for BHCs {fcoi = 1.7). 

Given that the hard flux is diminished greatly in the off/soft state of GRS 1758—258, it is 



likely that the Comptonization of disk photons is not as strong an effect as in other BHCs and the 
color correction may not be required. If no color correction is applied, for reasonable values of the 
black hole mass in GRS 1758—258 the inner extent of the accretion disk is broadly consistent with 
the marginally stable circular orbit around a Schwarzschild black hole (see Figure 4. This finding 
may be inconsistent with ADAF models of other BHCs at relatively low implied values of m (in 
some sources, however. X-ray flux and spectral states may be decoupled; see, e. g., Homan et al. 
2001; Wijnands & Miller 2001). If the color correction is accurate, the inner disk extent is slightly 
recessed and consistent with ADAF models of Cygnus X-1 in the low/hard state (rj„ > 20 Rg, Esin 
et al. 1998). However, the spectrum of Cygnus X-1 at lower implied values of m is more consistent 
with an inner ADAF volume than the spectrum of GRS 1758—258. A simple ADAF interpretation 
might still be valid for this source if the soft component could be attributed to an extended cool 
volume, but our results indicate that this is very unlikely. A model in which the inner advection 
region is peculiarly radiatively-inefficient might describe the off/soft state of GRS 1758—258. 

Smith, Heindl, and Swank (Smith et al. 2001d) have discussed the long-term spectral 
variability of BHCs IE 1740.7-2942, GRS 1758-258, GX 339-4, Cygnus X-1, and Cygnus X-3 
in terms of two independent accretion flows. As SHMS note, the off/soft state is adequately 
described by this model. The possible cooling of the accretion disk which we measure supports 
this picture. It is interesting to note that Chen, Gehrels, and Leventhal (1994) have proposed 
that IE 1740.7—2942 and GRS 1758—258 may accrete both through an accretion disk and via 
Bondi-Hoyle accretion from the relatively high-density ISM in the Galactic center. Observationally, 
the two-flow and disk-plus-Bondi-Hoyle accretion models may be very similar. 
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Table 1: Models and Results 









Multicolor Disk Blackbody Fits 






Model 
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Nmcd 


PL index 


NPL 
(10-3) 
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Blackbody Fits 
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Power-law Fits 
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(1022 cm-2 






PL index 


NpL 
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6 


2 q+o-i 






c 9+0.2 
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Note. — Results of fitting the background-subtracted, rebinned first-order RGSl and RGS2 spectra jointly with 
standard models, in the 0.6-2.3 keV bandpass (90% confidence errors). The power-law indices in models 1-5 are 
fixed to the values obtained by Smith et al. (2001) using RXTE within the off/soft state as this component cannot 
be constrained in the limited energy range of the RGS. The power- law component is not required in fits 1-5 and the 
quoted normalizations and fluxes are upper-limits (90% conf.). The neutral hydrogen column density, Nh, is allowed 
to float in models 1, 3, and 6. In models 2 and 5, Nh is fixed to the values reported by Smith et al. (2001) to 
allow for more direct comparison to those results. Similarly, Nh is fixed to the value reported by Goldwurm et al. 
(2001) in fits to XMM-Newton EPIC-MOS data obtained in September, 2000. Model 6 is a simple power-law; the 
obtained x^ statistic is slightly worse than the models dominated by thermal components, but the measured index 
is very different than those reported in other BHCs (see, e. g., Tanaka & Lewin 1995). In fits 1-5, the power-law 
contribution to the 0.6-2.3 keV luminosity is <2% of the total; luminosities are calculated assuming a distance of 
8.5 kpc. Luminosities quoted above are "unabsorbed" luminosities. 
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ulticolor Disk Blackbody Fit 
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Fig. 1. — Results of fitting to RGSl and RGS2 jointly with a model consisting of multicolor 
disk blackbody and power-law components. The flux (in units of normalized counts/cm^ /s) and 
data/model ratio are shown; at top: RGSl, at bottom: RGS2. Fits with a model consisting of 
simple blackbody and power-law components, and of only a power-law component (see Table 1) 
yield similar results in terms of the fit and data/model ratio. The large gap seen in RGSl (0.9- 
1.2 keV) is due to a failed CCD in that region of the dispersed spectrum. Points where the model 
goes through flux bins consistent with zero (in RGSl: near 1.2 keV, and 2.0-2.1 keV; in RGS2: 
between 1.0-1.3 keV, and 2.0-2.1 keV) are an artifact of our rebinning (requirement of 20 counts 
per channel) in regions of low or rapidly changing effective area. 
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Fig. 2. — Upper limits (90% confidence) on the strength of narrow, single-bin emission or absorption 
features in the RGS spectrum of GRS 1758—258. RGSl is shown in red, and RGS2 is shown in 
blue. The limits have been rebinned by a factor of 10 for visual clarity. These limits are based 
on our best fit to the continuum with the MCD multicolor disk blackbody model; results from a 
simple black body continuum model are very similar. 
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Fig. 3. — The combined RGS-1 and RGS-2 first-order lightcurve (time bins are 10 seconds, energy 
range: 0.3-2.5 keV). Above, the hghtcurve for the fuh observation; below, a randomly selected 
segment. Plotted through the lightcurve is the mean count rate (solid red line), and la deviations 
(individual error bars are la). Standard Shakura-Sunyaev (1973) accretion disks should vary on 
a viscous timescale (weeks). This lightcurve is extremely steady, supporting our interpretation of 
the low-flux state of GRS 1758—258 as being dominated by an accretion disk. 
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Fig. 4. — Inner accretion disk radii derived via the multicolor disk blackbody model, assuming an 
intermediate inclination (^inci = 45°) and a distance of 8.5 kpc. Fits with model 1 (Nh free) yield 
inner radii of 77 km and 220 km, for f^oi = 1-0 and fcoi = 1-7, respectively. Fits with model 2 
(Nh fixed) yield inner radii of 14 km and 39 km, for fcoi = 1-0 and fcoi = 1-7, respectively. The 
marginally stable circular orbit around a Schwarzschild black hole is 6 Rg; this radius is indicated 
with a dotted line. 



